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Abstract 

The D4-D8 brane system of Sakai-Sugimoto model at high quark density is studied in 
the weak couphng regime. We show that the color superconducting phase (for N(, 3) 
or the chiral density wave (for Nc oo) disappears at very large chemical potential, 
or equivalently at very large compactified dimension that the model possesses. We also 
comment on the prospects in the strong coupling regime along with the QCD phase 
diagram. 
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1 Introduction 



Decades have passed since quantum chromodynamics was accepted as the theory of strong 
interaction. While the successes of this theory are impressive, there still remain numerous 
unresolved problems. The most prominent of those problems is the confinement. At low 
energy, we observe hadrons instead of the QCD fundamental degrees of freedom and we still 
do not quite understand the mechanism of this phenomenon. The lack of the understanding 
is mainly due to the strongly coupled nature of the theory at low energy and to the fact that 
we do not have very good analytical control over the field theory in such coupling regime. 
The best hope, therefore, is that the numerical study could provide important insights into 
QCD at strong coupling. 

Confinement is one of the phases of QCD and there are other phases in different re- 
gions of the QCD parameter space. One such example is the deconfined phase at high 
temperature. The numerical method, in fact, is proven to be powerful in the investigation 
of the confinement /deconfinement phase transition, and has provided the estimates of the 
transition temperature and the order of the phase transition. (See Reference [1] and the 
papers cited therein.) There is, however, a large region of the parameter space in which 
even the numerical study has not been very successful. Namely, the theory at finite quark 
(baryon) density. There is a problem in carrying out lattice simulations with the nonzero 
chemical potentials that are conjugate to the density and this is known as the sign problem 
(see e.g. Reference |2]). While the effort to overcome the sign problem vigorously continues, 
there have been some developments in analytic approaches. These include the analysis at 
asymptotically large chemical potential and also the use of NJL type models. The former 
takes the advantage that the QCD coupling is expected to be weak at very high density and 
perturbative computations from the QCD Lagrangian itself are possible. The latter models 
are pure fermionic and quartic couplings, which reproduce properties of QCD interactions 
in a certain degree, are introduced^ 

What has been emerging from these analyses is the very rich phase structure of QCD 
in the parameter space of the temperature and the chemical potential. (See, for example. 
Reference [8] for the phase diagram.) In general, a cold and highly dense quark matter 
is expected to become a color superconductor [9]. As mentioned before, when the quark 
chemical potential fi is very large, the coupling g{fi) is small and the excitations near the 
Fermi surface of the quarks, particles and holes, are nearly free, and this naive ground 
state at high density is known as the Fermi liquid. However, the pairs of the particles at 
the antipodal points of the Fermi sphere are all degenerate and it costs no free energy to 
form such a pair. Then, if there is an attractive force between the particles (or holes), 
the pairing actually reduces the free energy of the system, leading to the instability of the 
naive ground state against the formation of the pair (Cooper pair). This is known as the 
Bardeen-Cooper-Schrieffer (BCS) instability [10] and in their original work for the electron 
gas in a solid, the attractive force was provided by the phonon exchange. For the case of 
high density QCD, there is an attractive force in a color channel and it leads to the similar 
BCS instability. 

In particular, if we have three massless flavors in the theory, a very interesting form of 

^ A good modern review on perturbative high density QCD is Reference [T]. The NJL-model was intro- 
duced by Nambu and Jona-Lasinio O U as the model that exhibits the chiral symmetry breaking and light 
hadronic spectrum. A good review on the NJL model is Reference [5]. The model is applied to QCD at 
finite density and reviewed in References [5] [7] . 
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the condensate has been found to form [lljrl This is the (scalar) diquark condensate of the 
form 

{QhQi,) = -iQhq'Rj) = ^CFLiStS', - 6^5t) , (1) 

where the superscripts a,b are the color indices, the subscripts i,j are for the flavors, the 
subscripts L, R indicates the chirality of the quarks and Aqfl is the size of the condensate 
(the gap). As one can see, the Kronecker deltas relate the flavor and color symmetries 
and those are not separately preserved by the condensate. The residual symmetry is the 
simultaneous flavor and (global) color rotations and for this reason, this phenomenon is 
called the color-flavor locking (CFL). The full symmetry breaking pattern of CFL is 

SU{3)coior X SU{3)l X SU{3)r x U{1)v x U{1)em ^ SU{3Uior+L+B X X U{1)q , (2) 

where S'C/(3)coior+L+H is the global diagonal subgroup of the original color and flavor sym- 
metries, Z2 is the subgroup of U{l)v that changes the sign of all the quarks and U{1)q is 
known as the "modified electromagnetism" whose gauge boson is a linear combination of 
the original photon and one of the gluons [ll^ [1] . We observe that since L- and i?-flavor 
symmetries both lock to the color, the chiral symmetry is broken through the color factor. 
Even though the mechanism of the chiral symmetry breaking is very unusual, the corre- 
sponding chiral Lagrangian can be built [12]. This novel phase of QCD created a renewed 
and wide interest in the QCD phase structure and many generalizations and modifications 
have been explored. For example, when one considers finite quark masses, there are other 
possible forms of the condensate and those can be energetically favored for some regions of 
the /i-T parameter space, resulting in the complicated structure of the phase diagram. The 
interested reader can pursuit the subject in the review papers cited above. 

Rather than continuing to overview the QCD phase structure, we would now like to 
turn to the high density behavior in the 't Hooft limit with weak coupling. (This is not 
QCD, which has Nc = 3, but has a potential relation to the holographic theories.) In 
this case, we do not expect the color superconductor to be the correct ground state of 
the cold QCD. Heuristically, this is because the Cooper pair is not a color singlet and not 
expected to survive the limit. The possibility of a color singlet condensate of particle and 
hole (not anti-particle) has been investigated by Deryagin, Grigoriev and Rubakov (DGR) 
in Reference |13j . When the particle and hole at the antipodal points of the Fermi sphere 
form a pair, the condensate is not homogeneous nor isotropic but is a standing wave in a 
certain direction; 

(g-L(x)(?R(y)) = e'P^<^+y^fix - y) , (3) 

where \pf\ = M ^^^d / is a function that describes the amplitude of the standing wave. The 
condensate is called the chiral density wave (xDW) and it breaks the chiral symmetry but 
not the gauge symmetry (in the limit a; — > y). This condensate is kinematically less favored 
than the Cooper pairing at A'^c ~ 3, but it has been shown that the ground state instability 
due to the formation of xDW (DGR instability) dominates over the BCS-type instability 
in the large Nc hmit [TSlITi] . 

The aim of this paper is to examine the high density behavior of the Sakai-Sugimoto 
model |il5t I16j . As we will describe in Section [21 this is a model in Type IIA string theory 
with a certain brane configuration. What makes this model interesting is that the low 

^ This model with Nf = 3 massless quarks is an approximate QCD where the masses of u, d, s are set to 
zero and those of c, b, t are taken to infinity. 
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energy spectrum is similar to that of QCD, especially, there are fundamental quarks with 
the left- and right-chiral symmetries. In the strong coupling gravity background analysis, 
Sakai and Sugimoto have shown that the chiral symmetry is broken and they were able 
to compute the hadron spectrum of the theory. Moreover, they have constructed the low 
energy effective action of the theory with the Skyrme term whose soliton excitations can be 
considered as baryons. Those promising successes make this model an interesting candidate 
for holographic QCD. Therefore, we are naturally motivated to examine the model at high 
density and ultimately, we hope the model to give insight into the structure of the QCD 
phase diagram in all values of the parameters, especially in the medium density region where 
the theory is strongly coupled. 

Though we do not get this far in this work, we show in Section [2] that the weak coupling 
regime of the theory at high density and zero temperature behaves very similar to QCD 
with a few differences. One of the differences is that when the chemical potential is very 
large with respect to the compactification size of the model, the BCS and DGR instabilities 
at A'^c ~ 3 and Nc ^ oo, respectively, are absent and the ground state is the Fermi liquid. In 
Section [Sj we comment on the known finite density analysis of the model at strong coupling, 
contrasting to the QCD expectations. In the strong coupling analysis, the possibilities of 
the superconductivity and x^W have not been addressed and we discuss some prospects in 
this direction. 

2 Weak Coupling Field Theory Analysis 

In this section, we discuss the Sakai-Sugimoto brane system at high density with the weak 
Yang- Mills coupling of the world- volume theory. Though the results are relatively straight- 
forward, the computations are somewhat involved. We therefore split the discussion into 
the qualitative and quantitative parts. In the first part, we qualitatively explain the high 
density behavior of the system, then in the second part, we carry out the computations and 
confirm the qualitative expectations. 

2.1 Qualitative Discussion 
2.1.1 BCS and DGR Instabilities 

Let us briefly review the BCS and DGR instabilities of high density QCD in a way that 
would provide the conceptual background for the quantitative calculations H For simplicity, 
we set all current quark masses to zero, which is a good approximation when the quark 
chemical potential is much larger than the mass of the heaviest quark. In the presence of 
the quark chemical potential, ^u, at zero temperature, we have a well-defined Fermi sphere 
of radius /i. As a convention, we take the excitations near the Fermi surface be particles as 
opposed to anti-particles. When is very large, the anti-particles are buried deep in the 
Dirac sea and will not play a role in the following discussion. 

On the Fermi surface, the free energy of the states are zero (more precisely, at the mini- 
mum) and it costs no free energy to change momentum along the Fermi surface. Therefore, 
the energy scales only in the radial direction of the sphere and we consider the renormaliza- 
tion group flow as we scale the energy down toward the Fermi surface. The relevant degrees 
of freedom are the particles and holes near the Fermi surface and we are interested in the 

For a modern exposition of the BCS instability, see Polchinski's TASI lecture notes [17j . 
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effective theory that describes the dynamics of those excitations. Let us specify a point on 
the Fermi surface by the momentum pp with the magnitude \pf\ = and decompose a 
four-momentum near this point as 

= {E, pF + l\\ + h) , (4) 

where is parallel to pp and l± is perpendicular to it. As stated before, only scales with 
energy and l± may be trivially integrated along the Fermi surface in a diagram computation. 
We, therefore, have the 1 + 1-dimensional effective theory that describes the dynamics of 
the particles and holes|l| It is important to notice that the kinematics is restricted because 
the dynamics must take place near the Fermi surface. When we consider a particle-particle 
or hole-hole scattering, it is clear that the scattering must be near back-to-back, that is, 
the scattering pairs must be at antipodal points of the Fermi sphere. In the back-to-back 
scattering, the scattering angle may be arbitrary without spoiling the kinematic restriction 
and hence the phase space of this scattering is all over the Fermi surface. 

Now in a two-dimensional theory, irrelevant operators of four dimensions may become 
relevant or marginal. In particular, a four-fermion interaction is marginal in two dimen- 
sions. If the interaction is attractive, the quartic coupling grows as we scale the energy 
down toward the Fermi surface and it eventually hits the Landau pole. This implies that 
the perturbation theory breaks down at the infrared scale around the pole. For the scat- 
tering of particle or hole pairs, this indicates that the naive ground state of the weakly 
interacting particles and holes near the Fermi surface, the Fermi liquid, is unstable against 
the formation of Cooper pairs. This is the BCS instability, and the new ground state has a 
gap due to the formation of the condensate whose size is roughly the location of the Landau 
pole. Technically, as we will see in the next subsection, the instability is closely related 
to the infrared divergence that appear in the perturbation theory and the gap properly 
provides the infrared cutoff. In weakly coupled QCD, i.e., QCD at very high density, the 
leading order contribution to the interaction is given by the one-gluon exchange and it is 
attractive in the antisymmetric 3-channel, resulting in the color superconductivity. 

We now turn to the DGR instability. This is associated with the scattering of the particle 
and hole located at the antipodal points of the Fermi sphere. In this case, the scattering 
is not back-to-back, but near forward. Unlike the back-to-back case, the scattering angle 
cannot be too large to stay near the Fermi surface and the phase space of the forward 
scattering is limited to a very tiny patch on the Fermi surface. The difference in the phase 
spaces for the BCS and DGR cases will be important in the quantitative computations. 

Now the particle-hole forward scattering is similar to the Bhabha scattering whose am- 
plitude has the forward enhancement. Thus also in our case, we can expect an infrared 
divergence (the DGR instability) as the exchange gluon becomes very soft, leading to a 
condensate of the particle-hole pair (xDW). This, however, does not happen in weakly cou- 
pled QCD. The reason is that the finite density screening effect completely overwhelms such 
a condensate; it is the screening that provides the infrared cutoff and not the formation of 
a condensate. 

The situation is different in the large Nc limit (with small 't Hooft coupling). First, 
note that the Cooper pair of the BCS instability is not color singlet while the x^W of 
DGR is. Therefore, the BCS instability is 1/Nc suppressed and DGR is not. Secondly, 

^ A rigorous derivation of the high density effective theory of QCD was carried out by Hong [18] . 
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since only the quarks are charged under the C/(l)v-symmetry, the finite density screening 
effect is provided by the quark loops, such as the one shown in Figure [H and the gluon 
loops do not contribute because the gluon propagator is independent of the quark chemical 
potential. This implies that the finite density screening is suppressed as the number of the 




Figure 1: One loop diagram that contributes to the finite density screening. 

color is taken to a large value. In fact, the DGR instability was discovered by disregarding 
the the screening effect and the authors of Reference |13j noted that the x^W can form 
and dominate over the Cooper pairing at least in the large Nc limit. Later, Shuster and 
Son showed that the DGR instability may occur if Nc ^ lOOOA'^^, where Nj is the number 
of the flavor [T^ . 

This is an example where the large Nc limit of QCD yields qualitatively different prop- 
erties. In weakly coupled QCD, xDW is not a relevant phenomenon. Though it could 
possibly compete with the Cooper pair at strong coupling, so far the situation is unclear 

m- 

2.1.2 Sakai-Sugimoto Brane System at High Density 

We now consider whether the Sakai-Sugimoto model at high density exhibits similar prop- 
erties as discussed above. For this purpose, we assume that the model is in the regime with 
low energy and weak Yang-Mills coupling so that we can use the perturbative world-volume 
field theory arguments. 

The model is Type II A string theory with D4-, D8- and D8-branes. The configuration of 
the branes is shown in Figure [2j The X4-direction is compactified to the circle of circumfer- 
ence L and the D8D8-branes are placed at the antipodal points of the circle. To discuss the 
low energy spectrum of the model, we first consider only the compactified Nc D4-branes. In 
Reference [20], Witten suggested that if we impose the anti-periodic boundary condition to 
the adjoint world-volume fermions, the low energy world-volume theory has the spectrum 
of 3 + 1-dimensional pure Yang-Mills theory. This is because the fermions get tree level 
mass of order 1/L and the scalars, including the compactified component of gauge field, A4, 
get the one-loop mass of order g/L, where g := g^l\fL and g^ := {2'K)'^gsls-, with gs string 
coupling and Ig string length scale, and we assume g <C lH 

Now, Sakai and Sugimoto insert the A'^^ D8 and Nj D8 branes as shown. As explained 
in Reference ^2]j|, there are massless fermions at the 3 + 1-dimensional intersection of D4- 
and D8-branes. These are the lowest states of the strings stretching from D8 to D4. Since 
the world- volume U{Nf)L gauge symmetry of D8-branes acts as the flavor symmetry, these 
massless fermions are fundamental "quarks". Similar massless fermions are also present 

^ Unlike four-dimensional case, other components of the gauge field do not acquire mass of order g^/L, 
in all orders of the perturbation theory. 
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Figure 2: Sakai-Sugimoto brane configuration. The a;4-direction is compactified with period L. The D4 
world-volume fermions have the anti-periodic boundary condition. We locate the Nf D8-branes and Nf 
D8-branes a,t X4 = = L and X4, = L/2, respectively. 

at the intersection of D4- and D8-branes. However, the GSO projection projects out op- 
posite chirahties to those fermions at different intersections. Therefore, we call the mass- 
less fermions at D4D8 and D4D8 intersections as "left-handed (ql) and right-handed (g/j) 
quarks" , respectively. 

This theory at low energy, therefore, has U{Nf)L x U{Nf)ji flavor symmetry and the 
theory appears to be very similar to the massless QCD, if Nc and Nf are appropriately 
chosen. The difference, of course, is that the fermions with different chiralities are separately 
located in the X4-direction and the gluons propagate in five dimensions, including the X4- 
direction. As argued by Antonyan et al. [22j, this QCD-like dynamics of the quarks at the 
intersections does not change even if the period L is large; the shift symmetry of the D4 
adjoint scalars and the five dimensional gauge symmetry allow the scalars, including A^, to 
couple to the fundamental fermions only through derivative interactins which is suppressed 
by the string scale. 

In the weak coupling analysis that we carry out in this section, we assume that the 
period L is much larger than the string length scale Is so that the tachyon becomes heavy 
and decoupleslfl Also as we have already explained, we place the D8 and D8 branes at the 
antipodal points of the compactified circle. 

In examining this model, we need to decide on how the quarks interact through the 
exchange of the gluons that propagate in the a;4-direction. One possibility is the non-local 
interaction. This scenario takes only into account of the zero- mode of the discrete mo- 
mentum in the X4-direction. In this case, the theory becomes completely insensitive to 
the existence of the fifth dimension and behaves in the same way as the four dimensional 

^ When the D4-branes of the system are replaced with their eflFective geometry (which is not the case in our 
discussion) , the proper distance of the D8s becomes less than the string scale in the region sufficiently near 
the horizon and the statement made here is no longer valid in the analysis with the background geometry. 
Refs.|23l 1241 125) show that the tachyon indeed condenses and it is responsible for the chiral symmetry 
breaking of the model at strong coupling. 
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QCD-like theory. We find this non-local scenario less appealing, especially when the com- 
pactification scale L is large. We thus take the second alternative where the D8-branes are 
treated as sources (or stiff walls) and allow the exchanged gluons to carry arbitrarily high 
momenta in the a;4-direction0 Such an assumption is reasonable because the D8 branes are 
infinitely heavier than the D4s and consistent with the fact that we are treating the D8s 
as the flavor branes, that is, we are neglecting their fluctuations. The momentum is not 
conserved in the j;4-direction but this is natural in that the translation symmetry is broken 
in this direction. 

We now explain how we introduce the quark chemical potential. We have the global 
symmetry SU {Nf)L x SU {Nf)ji x U{l)v- The chemical potential that we are interested in is 
conjugate to the U{l)v charge. The standard way to introduce a chemical potential in a field 
theory is to treat it as the constant background "gauged" field of a U{1) global symmetry, 
with all the components being zero except the time component. In this way, the chemical 
potential modifies the time component of the covariant derivatives in the Lagrangian of the 
theory. Thus in our case, the simplest way to introduce the quark chemical potential is to 
turn on the Aq constant background gauge fields of U{1) G U{Nf)L and U{1) G U{Nf)fi 
world-volume gauge symmetries and tune them to an equal value. Actually, the background 
fields may not be constant all over the D8-branes and the only requirement is to have the 
constant value at the intersections with the D4-branes. Thus, for example, we may turn on 
the field that depends on the radial direction in the 5,6,7,8 and 9 directions which would 
correspond to the nonzero electric field in the world-volume. 

What we will find in the quantitative analysis of the next section is rather intuitive. 
As we have mentioned, if we consider the effect only of the zero-mode momentum in the 
X4-direction, the theory reduces to the QCD-like theory. Thus when the compactification 
scale 1/L is very large compared to the energy scale of the interest, in this case it is the value 
of the chemical potential /x, we expect to have the BCS and DGR instabilities at A'c ~ 3 
and at Nc — > oo, respectively, just as described before. Now as the compactification scale 
1/L gets smaller, the infrared effect in the X4-direction becomes comparable to the one that 
leads to the regular BCS or DGR instability. As a consequence, the size of the condensate 
grows and eventually becomes too large to maintain the dynamics near the Fermi surface. 
Therefore, when the scale l/L is small with respect to fi, there is no BCS or DGR type of 
instability and the ground state of the theory is described by the Fermi liquid. 

In this qualitative discussion, it is not clear at what scale this crossover occurs. The 
computations of the next subsection show that at Nc ^ 3 and fiL > 1/g, no BCS-type 
instability is present and at A'^c — ^ oo and fiL > e^l^ j^fX with A := g^Nc, no DGR-type 
instability happens. Notice that the DGR instability persists to exponentially larger value 
of compared to the BCS case. This is because the phase space of the particle-hole 
scattering is very small, in fact it is exponentially small, and the discrete momentum in the 
compactified direction must become as fine as this scale to open up the extra dimension. 
The situation explained here is schematically summarized in Figure [3j 

The gluon momenta actually have to be cut off below the string scale to avoid the derivative interactions 
between the D4 adjoint scalars and the fundamental fermions. However, as we will see shortly, the high 
momenta increasingly suppress the gluon propagator and their contributions become negligible at sufficiently 
high scale. Therefore we can approximately take the a;4-momentum to infinity. 
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XDW 



BCS 



Fermi Liquid 



Figure 3: A schematic phase diagram of the theory at weak coupling, high density and zero temperature. 
Being schematic, the straight lines may not be straight nor sharp transitions in reality. 

2.2 Quantitative Discussion 

We now demonstrate quantitatively what has been discussed in the previous subsection. 
We carry out the renormahzation group and Dyson-Schwinger analyses. The former is 
more intuitive in accordance with the qualitative discussion and shows the existence of the 
instabilities. But this method does not provide the size of the gap and this is augmented 
by solving the Dyson-Schwinger equations. 

We adopt the conventions of Wess and Bagger [26], except the definition of the Dirac 
spinor; 



We mainly work in the chiral basis. As a convention, the undotted and dotted spinors live 
on the D8 and D8 branes, respectively. 

Our central focus of this subsection is to show the existence of the instabilities and to 
obtain the size of the gap. We are less interested in the exact color-flavor structure of the 
condensate, so in what follows, we simplify the analysis by suppressing the flavor structure. 
This is similar to Nf = 2 case where the Pauli principle requires the simpler quark pairing. 

2.2.1 Renormalization Group Equations 

Our first analysis is macroscopic in a sense that we introduce an effective one point four- 
fermion coupling. Then we observe how the effective coupling evolves as we scale the energy 
of the system down to the Fermi surface. This idea was first carried out in high density 
QCD by Evans et al. in Reference |27j . 

Since we are dealing with the weak coupling at high density, the quark interaction can 
be approximated by a single gluon exchange. We then model the four-fermion interaction 
by replacing the one- gluon exchange to a point, as shown in Figure HI Because the chemical 
potential breaks the 3-1-1 world- volume symmetry down to 0(3), we separately handle the 
couplings, and G-', as in 




(5) 



(6) 
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Figure 4: Replacing the one-gluon exchange to an effective one point interaction. 

For the one-gluon exchange, we can further write 

G\D) = -G\D):=-glX{D)F . (7) 

Notice that we have included the minus sign from the signature in the definition of . The 
constant is the five dimensional Yang-Mills coupling as before, F is the form factor that 
arises from the gluon propagator and 

X{D):=\{C{D)-2C{U)] , (8) 

with C{D) being the Casimir operator of SU{Nc) in the representation D and □ being the 
defining representation. 

We consider three color channels; symmetric (symm), antisymmetric (asymm) and sin- 
glet (•). For those cases, we have 

X(.yn,n,) = ^, X(^yu.n.) = . = , (9) 

where we have adopted the normalization tr(T'^T'^) = (1/2)5"'^ for the Nc x matrices 
{T" : a = 1, . . . — 1} of the defining representation. Notice the large Nc behavior of 
those factors. The singlet channel is larger than the other channels by a factor of Nc in 
the absolute value. This can be easily understood in the double line notation of the single 
gluon exchange diagrams, as shown in Figure [5l 




Figure 5: The double line notations of the one-gluon exchange. The top and bottom sets represent the 
symmetric (or antisymmetric) and singlet channels, respectively. If we fix the colors of the incoming quarks, 
the symmetric channel has the fixed colors for the scattered quarks, while the singlet channel has Nc choices. 

Let US now consider the form factor, F. The gluon propagator in the Feynman gauge has 
the form l/(p^+P4) with the discrete momentum, p4 = 2im/L, in the compactified direction 
and this propagator suffers from an infrared divergence. To cure this problem, we introduce 
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an infrared cutoff m. We simplify the situation by assuming that m is the same for the time 
and spatial components of the gluonslfl Possible one-gluon interactions are illustrated in 
Figure m In both cases shown in the figure, the discrete momentum is Inn/L because this 

0,1 




Figure 6: One-gluon exchange diagram. The gluons are propagating in the compactified 2:4 -direct ion. The 
D8D8-branes are treated as sources or stiff walls. 

is determined by the period of the compactification and not by the distance between the 
branes. Therefore, the form factor, that arises from the gluon propagator, is independent 
of the D8 brane distance and is the same for the both interactions in Figure [6] (including 
the case with the gluons making many rounds in the compactified circle). We follow Evans 
et al. [23 to obtain the form factor, namely, we take average of the gluon propagator over 
the scattering phase space. For the dynamics very near the Fermi surface, the energy and 
the component of the momentum in the radial direction of the Fermi sphere are almost 
zero. Hence apart from the extra X4-direction, the phase space is two dimensional along the 
Fermi surface. Let pi and p2 be the incoming and outgoing momenta of a scattering quark. 
We then have the gluon momentum 

= {pi-P2? -2i?{l-cose) , (10) 

where is the scattering angle. 

For the back-to-back scattering, the angle ranges from to vr and the phase space is all 
over the Fermi surface. Therefore, the form factor in this case is 

1 ^ /■ d?p 



n=— 00 

27r 



NL ^ J p'^ + (27rn/L)2 + m 



2 



NL 



In 



sinh (^mL^Jl + A^j? /m?^ / sinh ^—mL 



(11) 



where we have defined the total phase space factor J\f := Airfi^. We have already carried 
out the sum over the p4-discrete momentum because the dependence is only in the gluon 



* At very high density where the coupling is weak, one can expect the finite density Debye screening of 
order provides the infrared cutoff. This is true for the time component but not for the spatial components. 
There is no static screening in the spatial (magnetic) components and the magnetic screening is dynamical 
due to the Landau damping. This was pointed out by Son [28] and he discovered that the dynamical 
screening effect leads to the qualitatively different form of the gap. We will take into account of Son's effect 
in the Dyson-Schwinger analysis. 
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propagator. As an effect, we have encoded all the information about the extra dimension 
into the form factor. For the forward scattering of a particle-hole pair, the angle should 
not take very large value so that the dynamics takes place near the Fermi surface. Thus 
we require, in this case, the angle range from to ^c/y, where the latter angle is much less 
than vr and limits the phase space to a little patch on the Fermi surface. The form factor 
then takes the form 



Fpw 



1 



ML 

27r 
JdL 



E 



In 



+ (27rn/L)2 + m? 
sinh ^-mL\/l + 2(1 — cos 0;7y)/x^/m^^ / sinh ^2"^^^ 



(12) 



where M := 27r/Li^(l— cos Ouv)- Notice that since Ouv ^r, we have MFpw < MFbb- This 
fact will lead to the dominance of the BCS-type instability over the DGR type for Nf. ~ 3. 
When the dimensionless parameter mL is small, the form factors behave logarithmically 
with respect to the parameter fi/m and they behave linearly when niL is large. This is 
because when the the compactified dimension is very small, the contributions from n 7^ 
is also small and the integral is effectively two dimensional, leading to the log behavior. 
When the compactification size is very large, the discrete momentum becomes finer and the 
integral essentially becomes three dimensional and the form factors behave linearly. 

Having modeled the four-fermion interaction, we now derive the renormalization group 
equations for the couplings. The diagram that drives the renormalization group flow is the 
fermion one-loop diagram and we consider the three cases shown in Figure [71 For Diagram 




R 




R 




(a) 



(b) 



(c) 



Fi gure 7'. The diagrams that drive the renormahzation group flow. The letters "Z/" and "i?" refer to left- 
and right-handed quarks, respectively. 



(a) in the back-to-back scattering, we can deduce from the expressions 
corresponds to either 

, or . 
When the both vertices correspond to G^^, the diagram yields 



that each vertex 
(13) 



05a -O-y I3\ 



PIS 



{px - n6x,of 



{-Px - l^hflY 



0«(5-0/37 



), (14) 



where we used the quark propagators shown in Appendix [Al Note that the momentum 
integral is four dimensional rather than five. The dynamics of the quarks are restricted 
to the four dimensional intersections of the branes and the information about the extra 
dimension has been encoded in the form factor. We decompose the momentum as in Equa- 
tion dH) but redefine pp to include l±. Then near the surface of the large Fermi sphere, the 
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vectors pp and /y are near parallel and we also have E^l\\ <C /x and pp ^ A*- Under these 
approximations together with the use of the 0(3)-invariance, the argument of the square 
bracket in Equation ()14p becomes 

One can integrate over E, either by the contour integral or by the Wick rotation, then |/||| 
is integrated from the scale Ajjv down to A/^. One can also simplify the cr-matrices (all the 
necessary formulas are given in Appendix B of Wess and Bagger [26]) and the expression 
(fT^ becomes 

^Ml)' - l^^''-^''-') t , (16) 

where we have defined t := ln(A//j/A{/v'). This parameter t has the range (— oo,0) and the 
lower limit corresponds to the Fermi surface. When one of the vertex of Diagram (a) is 
and the other is , similar procedure yields 



I. 



(GIlGD (-2cT0^^a0^^ + ^a^^^a^^^^ t , (17) 



167r3 

and when the vertices are both G^ , we get 



(Giif (^5a°'^'^cT°^T - H^i55-i77^ t . (18) 



167r3 

From those results, we obtain the renormalization group equations 



dGlp 


= ^ 1 


r_ 


dt 


167r3 1 




dGi, 


= \ 




dt 


167r3 1 


[3 



{Gl,f + 2Gl,Gi,-5iGi,f} 



^2 10 13 
These equations can be diagonalized to the following forms 



(19) 



dt 167r3 



{GIl - ^Gi^Y , (20) 



It --^(Gpp + G^l) . (21) 

One can carry out the same procedure for the back-to-back scattering of Diagram (b) in 
Figure [7] and obtain 

d(G^r p + 3G{ p) , , 

dt 

dt --^^^LR-Gpji) . (23) 

Above four renormalization group equations are obtained by Evans et al. [27] (but with 
different form factors). 
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The Diagram (c) of Figure [7] is similar to the case with (a) [and not with (b) because qr 
has an undotted spinor index just as ql]. But we should recall that we are interested in the 
forward scattering in this case. Therefore in the loop, the top left-handed quark propagator 
carries the momentum (E, pp + Z||) and the bottom right-handed one carries {E, —pp + 
in the directions of the arrows in the quark loop. Then the propagator part of the diagram 
corresponding to (fT5|) of the Diagram (a) is 

- i ("°-^°'' - r'""^'"") ^ • 

This structure is the same as the LL-case (a), except the overall sign. This sign is cancelled 
by the other one that comes from the difference in the direction of the quark line. Thus, 
apart from the integration range of the scattering angle 9, the cases (a) and (c) are the 
same. So we have for Diagram (c), 

■^^^^ = -T^Kn-^oi,f. (25) 

(20) 



Now for the generic form of the renormalization group equation 

dG{t) 



-KG(ty , (27) 
with some constant K, the Landau pole, if exists, is reached at 

Since the range of t is (— oo,0), the Landau pole exists only when KG{0) is positive and 
the larger the factor KG{0) is, the faster the pole is reached. We note that because of 
the constants M and M in the form factors (jlip and (jl2p . the Landau pole is independent 
of those. From Equations ([7]) and ([9]), and also with the fact that AiFpw < J^Fbb, 
we see that the instability is dominated by the BCS type in the LL, color antisymmetric 
channel for ^ 'i. Note that the color symmetric channel does not have instability 
because the interaction in this channel is repulsive. When Nc is sufficiently large, the LR- 
channel, whose coupling is proportional to the product of X{») and Fpy/, dominates over 
the other channels because the inequality M.Ffw < J^Fbb is compensated by the fact that 
|X(asymm)| < Thus the DGR-type instability dominates over BCS in this regime. 

In our method here, it is not possible to estimate the value of Nc at which the crossover 
from the BCS- to DGR-type instability occurs, because we crudely introduced the common 
infrared cutoff, m, to all the components of the gluon propagator and we do not have the 
actual value of Ouv See Shuster and Son [H] for the estimate of the value A'^c- 

Recall that the form factors grow linearly when the parameters mL and p/m are large 
and the effective four-fermion couplings become large accordingly. We thus expect that the 
whole analysis breaks down when those parameters are exceedingly large, and we need to 
resort to the microscopic analysis, i.e., the analysis with the fundamental interactions, to 
gain insight into the nature of the pathology. This is the subject of the next analysis. 
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2.2.2 Dyson- Schwinger Equations 

We now turn to the analysis based on the Dyson-Schwinger equations. The traditional form 
of the Dyson-Schwinger equation in the diagrammatic representation is shown in Figure [8j 
This method is less intuitive compared to the previous renormalization group analysis, but 




Figure 8: The diagrammatic representation of the Dyson-Schwinger equation for the diquark condensate. 
The triangle denotes the gap insertion. The square in the gluon propagator represents the screening effect. 

it lets us obtain the actual size of the gap. This analysis is microscopic which deals with 
the quarks rather than the quasi-particles of the effective theory and interactions are the 
QCD interactions rather than the effective ones. 

As usual, we adopt the Nambu-Gor'kov formalism (see, e.g., Reference [9J). In order to 
introduce the Nambu-Gor'kov basis, we define the charge conjugate Dirac spinors as '0'^ '■= 
C'ilF and := ijFC'^ where the charge conjugation matrix C is defined in Appendix Rl 
In the Weyl basis, these can be expressed as 

(29) 

Then we define the Nambu-Gor'kov basis as 

^ - « =:c ^I3\'^ ,T, ._ 1 a cf3 



G{p)- 



(32) 



^ := ^ [qLc., Qr'', qlp, Cj ' ^ ^"^^ '^L ^ Irp) ■ (30) 

The advantage of the Nambu-Gor'kov formalism is that we can naturally include the con- 
densates in the propagator of ^. For example, the diquark condensate in the s-wave (LL 
or RR condensate) is given as 

V^^C^^S^ = -iqiQL + iq^R , (31) 

thus the inverse propagator that contains this condensate can be written as 

/ {p^ - f^S^,o)cT'' iAnip) \ 

(Pu - 1^6^,0)^'' -iAL{p) 

iAL{p) {p,. + i^6^,o)(t' 

\ -iAnip) {Pu + f^S^fi)^'' J 

The A-matrices appearing in the inverse propagator are defined as 

AlMp) = A+{p)Pl,r+{p) + A^{p)Pl,r-{p) , (33) 

where A±{p) are the gaps and the quark (anti-quark) on-shell projectors, Pl,r±{p)i are 
defined in Appendix |Al The ones with the bar can be obtained by replacing P — > P in the 
above expression. (We assume the gaps to be real.) We can invert the matrix (j32p by using 
the formulas listed in Appendix lAI^I If we write 

G{P) = {^'' ^'') , (34) 



® Note that the projectors are not invertible, so one must use appropriate inversion formulas. 
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we then have 



A_P 



G 



21 



-"-(Ip1-m)2-a^ 







A+Pfl_ 





+ 



A-Pfl+ 
pg-(|pl+/^)"-A^_ 



(35) 



and other components will not be important in writing down the Dyson-Schwinger equa- 
tions. Notice that if the condensates, A-t, in the denominators vanish there are terms that 
diverge as the energy is scaled toward the Fermi surface, i.e., po ^ and \p\ ^. This is 
essentially due to the sign structure of the chemical potential in the matrix (j32p and to the 
fact that the A-matrices occupy off block-diagonal components. This infrared divergence is 
cured by the formation of the condensate and the A properly behaves as such a condensate. 

If this were the traditional four-dimensional set up, we could have introduced the p-wave 
diquark condensate (Li?-pair) in the Nambu-Gor'kov propagator. [Such condensate would 
have occupied the anti-diagonal slots of the matrix (j32p .] However, this is not allowed in 
our theory. As shown in Figure EJ since ql and qn separately live on the D8 and D8 branes. 



Allowed 



Not AUowd 





D8 D8 D8 




Fi gure 9* The left diagram represents the right-hand side of the Dyson-Schwinger equation (Figure [Sj for 
LL-diquark condensate. On the right panel is the similar diagram with I/i?-diquark condensate. This clearly 
is not making sense, because the condensate (the triangle) is separated and also is not being able to flip the 
helicity. 

the Dyson-Schwinger equation for such condensate cannot make sense. In the previous 
macroscopic renormalization group analysis, we encoded all the information about the extra 
dimension in the form factors and the left- and right-handed quarks effectively lived in the 
same four dimensional spacetime. However, in this microscopic Dyson-Schwinger analysis, 
we see that it is actually not possible for the -Li?-condensate to form. We emphasize that 
the condensate is not being energetically suppressed but simply not possible to form in the 
brane system of Sakai and Sugimoto0 

For the particle-hole {LR or LR) pair, we might naively introduce the condensate in 
the diagonal slots of the matrix in (j32p . However, this represents the introduction of the 
usual chiral condensate, that is, the pair of particle and anti-particle, and not the desired 
particle-hole pair. One can verify in this case that the infrared divergence near the Fermi 



Refs. |29l 130) discuss the condensation of a gauge invariant particle-antiparticle pair across the D8 and 
D8 branes at strong coupling. Although it is possible that this sort of condensation forms in our situation, 
we note that we are considering the condensation of a particle pair which is not gauge invariant and at weak 
coupling. Therefore, this possibility is not taken into account in our discussion here. 
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surface is not present. Thus, the Nambu-Gor'kov basis is not suitable for introducing the 
particle-hole chiral condensate. Instead, we propose to use the doubled basis. 



[iLa, qn , QLa, QR ) 



(36) 



and consider the inverse propagator of the form (j32p with the replacements, iAl — > S/,, 
—lAji Tjji, —iAl '^l and zA/j — > S^, where the S-matrices are similarly defined as for 
the A-matrices. The spinors of the second set have the chemical potentials in their kinetic 
terms with opposite sign from the first set. This effectively introduces the hole degrees of 
freedom. In this case, the propagator has the infrared divergence near the Fermi surface 
and the S-condensate provides the cutoff. Thus we properly have the interpretation that 
the condensate is the particle-hole pair near the Fermi surface. We note that since the 
condensate is formed out of the spinors with the dotted- or undotted-index pair, and not 
the mixed one, the condensate lives either on D8 or D8 branes and not across them. 

We must also consider the gluon propagator which is the other ingredient of the Dyson- 
Schwinger equation. Unlike previous macroscopic treatment, we properly take into account 
of the perturbatively computable screening effect. The most general form of the 0(3)- 
invariant gluon propagator is 



+ 



p^Ap) 



p2 + (27rn/L)2 + G,(p) p"^ + {2TTn/L)^ + Fs{p) 



(37) 



where Fg and Gg are the electric and magnetic screenings, respectively, and the projectors 
are defined as 



Pijip) = Vij 



PiPj 



pUp) = = pUp) , pUp) = V^u 



p2 



flU 



(38) 



We have dropped the term with the gauge fixing parameter in the propagator. This term 
has been verified not to contribute to the gap, the solution to the Dyson-Schwinger equation, 
at very large chemical potential |31j . 

As explained in Section 12.11 the finite density screening effect is given by the diagram 
shown in Figure [H at one-loop level. Since the quark loop stays on the D8 or D8 branes, 
there is no extra dimensional effect on the screening, and the gluons with nonzero momentum 
in the j;4-direction, such as the case shown in Figure [9l do not have the screening. We thus 
have the standard expressions [32] 



Fg{p) = 2ml^ 
Gs{p) = mlj^ 



1 _ ( El 

\p\ vIpI 



Qoix) 



1 



In 



l + x 



1 — X 



vr 



^n,0 
PO 



i-e{l-x^), ml 



1 

4^ 



(39) 



where 9 is the Heaviside function and 6nfl signifies that the screening is effective only for 
the gluons with n = 0. 
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We can now write down the Dyson-Schwinger equation. We start with the equation for 
the A-condensate. At high density (weak coupling), we can approximate the quark-gluon 
vertex with the bare ones and we have 



G{k)-'-Go{k) 



-1 



-ig^X{D) 



T^G{p)T''D^,{q,n) . 



(40) 



We have included the color factor X{D) as in Equations ([9]) and the representation D is 
either symmetric or antisymmetric depending on the channel that G is in0 Also we have 
defined Go as G without the condensates, q := k — p, and are 



:= 










^ 































/ 



(42) 



We can use (I32]l . (pjSl) and (j40l) to obtain the gap equation 

d'^p 



A±{k) =ig^X{D)S^ 



n 



tr 



(2vr)4 

A+(p)Pfi,(p) 
-{\p\-^,f-. 



+ 



A_(p)Pfl+(p) 



a''PL±{k) 



(43) 



where the trace is over the spinor indices and is taken to project out A± from the A^-matrix 
in Equation (j33p . After some algebra, one obtains 



A 



n 

A+(po 



d^p 



pI 



+ 



-^)2-A+(po)2 
A-(po) 



1 - {p-q){k ■ q) 



+ 



h + lP-k 



Pi 



q^ + {2TTn/Lf + Gs{q) q^ + {2TTn/ Lf + Fs{q) 
l + {p-q)ik-q) , \-\p-^k 



+ 



+ ^iY-A,{p^YW + {2TTn/Lf + Gs{q) q' + {2TTn/ Lf + F^iq) 



(44) 



where we have assumed that the gaps are functions only of fco or pQ. In deriving this 
equation, we have adopted the approximation, qo \q\ ^ /i, so that P^^, rj^^df^fldufl- 
The equation for A_(A;o) is the same except that the two terms in the round brackets are 
exchanged. Note that only the first term in the equation of A_|_(A;o) has the near-Fermi- 
surface (infrared) divergence that is being cured by the formation of the condensate. Thus 
to the first approximation in large /i, we can neglect the second term. Similar observation 
in the equation of A_(fco) results in the conclusion that this gap does not form at near the 
Fermi surface. 



The easiest way to see how this color factor comes in is to note that 

^T^fcTc"^ = ^X{symm){5ab5cd + SadSbc) + ^X(sisymm) {5 abScd — SadStc) 



(41) 
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Now as before, we set p = Pf + and approximate \pf\ ^ pf ■ l\\ ~ 1, |pl ~ /i + l\\ 
and q'^ ~ 2^'^{\—p-k). Then the integration measure takes the form fj^dpodl^^d cos 6d(j) with 
cos 9 := p ■ k. Since the integral is dominated by the region ~ 0, we further approximate 
that p ■ k ^ 1 and p-q^O^k-qin the numerators of the integral, but not in the 
denominators. We Wick rotate po ipo and integrate over /y and (p. The gap equation 
now takes the form 



A+(A:o) ~ / dpodcosei^- 



cos e + (l/2){27rn/(^L)}2 + Gs{q)/i2fi2) 

, 1 \ A+(po) 

+ 1 _ cose + (l/2){2vrn/(^L)}2 + FM/i^f^^)) VpI + A+bo)^ ' ^ ^ 

with the approximate form of the screenings 

F,{q)^2ml6n,o, ^ ^m?,^ 5„,o • (46) 

From this equation, it is clear that for the symmetric channel D = symm, we only have the 
trivial solution A+ = 0@ We thus consider the antisymmetric channel from now on. 

In Equation (j45p . the sum over n and integral over 9 can be carried out in a straight- 
forward manner and yields 

A.(^o) - / dp. In f — m ^ , (47) 



where we have defined 



A := 2^0^/2^^ Af^^/^(7-V{smh(/iL)/(/iL)}4 , (48) 

and the part, {• • • is the contribution from the sum over n 7^ 0. The factor |A;o — Po\ 
appearing in the logarithm comes from the Landau damping of Gs{q) as in ()46p and this 
effect was first discussed by Son [28j. We can follow Appendix B of Reference ^28j to solve 
this equation and obtain 




A+(fco) = Aosin [ ^/ _(_\ 7~ 1 ' ^i*'^ Aq = Aexp 





127r2 2Nc 



(49) 

Let us comment on this result. We first note that the gap vanishes in the 't Hooft 
limit. This is consistent with what we have concluded in the renormalization group anal- 
ysis. Now, when /xL ^ 1, we have smh{^L)/^L ~ 1, so the extra dimensional effect in 
A disappears and the resulting expression for the gap coincide with the QCD result for 
Nc = 3. (See, for example. Reference [1].) In the opposite limit where fiL — > 00, we have 
sinh(//L) exp(//L)/2, so the extra dimensional effect contributes heavily and the gap 
grows with the parameter ^L. As one can observe in Equation (j45p . this is because the 
infrared effect of the terms with n 7^ comes to be comparable to that of the n = term. 



If the flavor structure is included, this conclusion gets slightly more complicated. However, the fact 
that the antisymmetric channel dominates over the symmetric one does not change. See, for instance. 
Reference l33l. 
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From Equation ()49p . we see that at /xL ~ the gap is much larger than the size of the 
Fermi sphere itself and such a solution cannot be accepted, for the dynamics is no longer 
taking place near the Fermi surface. We therefore conclude that when fiL ~ 1/g, the gap 
does not form and the ground state simply is described by the Fermi liquid. 



We now turn to the S-condensate, xDW. The computation is almost identical to the 
previous case and we arrive at the equation similar to Equation (j45p with the replacements 
A+ — > and £)—>•. There are, however, a few differences. The most important one is 
the range of the integration parameter 6. This is restricted to the near infrared region, i.e., 
6^0, because this is a forward scattering and the exchanged gluon should not be harder 
than the size of the gap or the momentum carried by the propagator G{k). When the angle 
is set to small value, the /n component is about — cos 6), and the propagator carries the 



momentum approximately a/pq + . We thus require 



/x(l-cos0) < ^p^ + S^ , (50) 

and this inequality sets the upper limit on 6. Now because of this kinematic restriction, when 
Nc ~ 3, the A-condensate dominates over the S-condensate. Therefore, in the following, 
we consider the 't Hooft limit (with small A). In this limit, the screening Fg and Gg is 
1/iVc-suppressed, so we drop the screening terms from the gap equation. In the absence of 
the screening, the approximation, q'^ ~ 2//^(l — cos9), has infrared problem when n = 0. 
This means that Qq cannot be neglected in this case and we must use ^ 2/i^(l — cos 9) + 
1^0 —pol'^^nfi- We thus have the gap equation for the S-condensate 

S+(A;o) =-§^\^Yl J dpodcosO 

2 S+(po) 



I- cose + \ko - Pol V(2/x2)5„,o + (l/2){27rn/(^L)}2 ^p^ + S+(po)2 



(51) 



where the integration range of 9 is restricted as mentioned above. 

Let us first consider the case with n = 0. In this case, we can carry out the cos 9 integral 
with the restriction (1501) and obtain 



We can again solve this equation following Son [28j (also see [I9j). In this case, the cal- 
culation is slightly different from Reference [28], so it is shown in Appendix [Bj The result 
is 



E+(fco) = So cos [ \l ^ ^"jy ^ ^'^ I ' "^it^ So = 2^exp 



/47r2 Nc 



(53) 

This result agrees with Reference jl3j . 

When we include the terms with n 7^ 0, the logarithm term in Equation (j52p gets 
augmented as 
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where e := \Jp'q + If ^.L is not too large compared to e^/^, then the second term in 

the above expression is small and the result (j53p does not change. However, when /iL is so 
large that the second term yields dominant contribution ~ /iLe^/^, the integrand of the gap 
equation (j52p takes the from proportional to l/(po + S^)^/^, which is free of the infrared 
divergence even without the condensate S. This implies that the gap does not exist. To 
roughly estimate the value of fiL at which the crossover occurs, we approximate the first 
term in Equation ([M]) as ln(2/x/So) and the second term as /xI/(Eo//i)^/^. Then we see that 
the second term becomes important when 

> (/x/So)i/2 ln(/x/So) « e^'^/VX . (55) 

3 Comments on the Strong Coupling Regime 

We have addressed the possibilities of the color superconductivity and the chiral density 
waves in the Sakai-Sugimoto model at finite density and explicitly carried out the compu- 
tations in the weak coupling limit. As was stated in the introduction, the ultimate goal 
of this investigation is to obtain the quantitative behavior of the model in the weak and 
strong coupling regions, make qualitative comparison of the phase diagrams and gain insight 
into the QCD phase diagram which is still unsettled. We therefore comment on the strong 
coupling gravity background analysis of the model at finite density. 

As suggested by Sakai and Sugimoto [15], the strong coupling analysis is done by tak- 
ing the gravity background limit of the D4-branes while treating the D8-branes as probes, 
then the DBI-action of the probes are studied to obtain the spectrum of the low energy 
excitations at strong coupling. The generalization to the finite density has been discussed 
in References |341 [35} 136] . The phase diagram of the model in the space of temperature 
and chemical potential has been obtained by Horigome and Tanii [35] . Let us briefly review 
their results. In the previous section, we have placed the D8- and D8-branes at the an- 
tipodal points of the compactified circle. However, this is not necessary and in this gravity 
background analysis, the compactification radius is set to R and the distance between the 
branes to L with the range < L < irR. 

At finite temperature, in addition to the X4-direction, the Euclidean time direction, r, is 
also compactified and the period of the time circle is identified with the inverse temperature. 
Horigome and Tanii consider the three known phases, first discovered by Aharony et al. |37] 
at zero density, and the spacetime configurations for the phases are illustrated in Figure [101 
In the figure, the vertical [/-axis represents the radial direction in 5,6,7,8 and 9 directions. 
The change in the background geometry of the phases from (a) to (b,c) is interpreted as 
the confinement /deconfinement phase transition by Aharony et al. The thinner lines and 
curves in the diagrams represent the D8 or D8 branes and the change in the configuration 
from Diagram (b) to (c) is interpreted as the chiral symmetry restoration. 

The on-shell DBI actions of the D8-branes with nonzero chemical potential have been 
obtained by Horigome and Tanii for each phase. For the configurations with the smooth 
U-shaped D8-branes, that is, for the diagrams (a) and (b) of Figure [TOl they have found 
that the chemical potential must be constant along the radial [/-direction and the actions 
for those configurations are independent of the chemical potential. Only for the parallel D8- 
brane configuration of Diagram (c) has the non-trivial dependence on the chemical potential 
in its action. By comparing the on-shell actions at the various values of the temperature and 
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(a) confined ^SB 



(b) deconfined ^SB 



X4 



(c) deconfined 



Figure 10: The spacetime configurations of the three phases. The vertical axis U is the radial direction 
in 5,6,7,8 and 9 directions. The thinner lines represent the stacks of D8 and D8 branes. Diagram (a) is 
the low temperature confined phase. The chiral symmetry is broken in this phase. Diagram (b) is the high 
temperature deconfined phase and with xSB. Diagram (c) is also the high temperature deconfined phase 
but with chiral symmetry restored. 



the chemical potential, they determined the phase diagram which is schematically shown in 
Figure EIiI 



Deconfined 
Deconfined ^SB 




Confined ^SB 



(a) L>0.97R 



(b) L<0.97R 



Figure 11: The schematic phase diagram obtained by Horigome and Tanii. The temperature Tc denotes 
the confinement/deconfinement phase transition temperature. The inter-D8D8 distance L ~ 0.97/? is the 
critical value where the deconfined phase with xSB exists. 

The confinement/deconfinement phase transition line aX T = is determined by the 
D4-background geometry and the D4 difference action in those phases scales as N'^ |37j . 
This is expected to completely dominate the D8 probe actions, which scale as NcNf, so 
the confinement/deconfinement phase transition line is not affected by the value of the 
chemical potential^ Therefore, Horigome and Tanii assume that the phase below the 
confinement/deconfinement line is in the configuration (a) of Figure [TOl and compare the 
D8 actions of the configurations (b) and (c) in the deconfined D4-background geometry 
to obtain the phase structure above the confinement/deconfinement phase transition line. 
This is why the xS/xSB phase transition line of the right panel in Figure [11] is terminated 
at r = Tc- When the inter-D8D8 distance L is larger than 0.97i?, it has been shown by 
Aharony et al. [37] that the phase represented in Diagram (b) of Figure [TOl does not exist. 
Therefore, the phase diagram becomes rather structureless, as shown in the left panel of 
Figure [TTl 



In Reference [35], the temperature and chemical potential are measured in different units. If measured 
in the same units, the chemical potential is asymptotically larger than the temperature by the factor of the 
't Hooft coupling A. 

In the original version of this e-print, we overlooked the dominance of the D4 action and compared only 
the D8 actions in the confined and deconfined phases, which was incorrect. 
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The analysis of Horigome and Tanii assumes only three possible phases and this is similar 
to the simplified picture of QCD where one assumes only hadronic and quark-gluon plasma 
phases. In the latter, it is expected that when the chemical potential is raised several times 
h.QCD (at low temperature), the hadrons start to overlap and the quarks are shared by 
many hadrons. This implies the change in the degrees of freedom and since only the two 
phases are assumed, there should be a phase transition. Stated differently, in this simplified 
picture of QCD, we expect that a phase transition to occur even at zero temperature. See 
Figure 1121^^1 Though the physical setups are similar, we see that the phase diagrams in 



T 





Qaurk Gkion Plasma 






Hadron \ 






Figure 12: Expected phase diagram of the simphfied QCD where only the hadronic and quark-gluon 
plasma phases are assumed to exist. 

Figures [TT] and [T2] are qualitatively different. We suspect that the discrepancy stems from 
the large Nc limit and the probe approximation in the Sakai-Sugimoto model. As we have 
explained, the confinement /deconfinement phase transition is determined by the D4 actions 
because the difference action scales with A'^ and dominates over the D8 actions (with the 
chemical potential). It is unlikely that this picture can change in the probe approximation. 
Therefore, in order to see modifications in the confinement /deconfinement phase transition 
line, one should take into account of the D8-brane back reaction to the geometry. 

It is also clear that we do not observe the color superconductivity or chiral density waves 
in the strong coupling analysis of Horigome and Tanii because these possibilities are simply 
not considered. To explore those exotic phases, one must come up with the corresponding 
stringy pictures of the branes and strings, and see if they are energetically preferred at any 
point in the /i-T parameter space. Such stringy pictures of quark matter are, so far, not 
clear to usJifl Nevertheless, we make some general remarks in this direction. First, we must 
remember that the relevant degrees of freedom at high density are particles and holes while 
the anti-particles are buried deep in the Dirac sea. Thus, if the U-shape configuration of 
the D8-branes describes the mesons, which are the pairs of particle and anti-particle, then 
we expect the configuration of the superconductor or chiral density waves to be different 
from the U-shape configuration because they are the pairs of particles and holes. Secondly, 
the Fermi sphere plays the essential role in high density QCD. Thus, the Fermi sphere must 

It is interesting to observe that the weakly coupled large Nc Af = 4 super- Yang-Mills theory on three 
sphere and Type TIB supergravity on AdSs x both have qualitatively the same phase structure as Figure [T2l 
[381 1391 140| . In these cases, however, the chemical potential is conjugate to the 17(1) subgroup of SU{4:) 
7?-symmetry. 

The stringy configurations of baryons have been suggested in Refs. [411 142) in which the nuclear matter 
phase has been discovered in the phase diagram. 



23 



be encoded in the holographic picture somehow. Also, the holography is in the 't Hooft 
limit under which we do not expect non-gauge invariant quantities to survive. Therefore, 
it is likely that we do not observe the color superconducting phase and the relevant phase 
probably is the chiral density waves in the holographic theories. 

Although Shuster and Son [13] settled that the color superconductivity dominates over 
the chiral density waves in the high density weakly coupled QCD (with Nc = 3), it still is 
not known if this observation persists at the medium density strong coupling region. The 
Sakai-Sugimoto model in the usual 't Hooft limit is not likely to be able to address this 
problem but it is interesting to examine this in the limit with Nf /Nc fixed. In this case, we 
can expect the effect of the color-flavor locking to be significant and also the competition 
between the D4 and D8 actions becomes non-trivial. 

Finally, as we have mentioned in the introduction, the rich structure of the QCD phase 
diagram is partly due to the quark masses. As the authors of Reference |21) have already 
noted, the inclusion of the quark masses involves the tachyon that comes from the string 
stretching between the D8 and D8 branes. This subject is being actively studied (for 
example, see Refs.[231 [Ml [3D]) and the resulting phase diagram is yet to be seen. 

//the finite density holographic model turns out to capture the aspects of QCD, it would 
be very interesting because we can explore the region of the QCD phase diagram where the 
perturbative nor the numerical analysis is available. 
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A Some Formulas 

We adopt the convention of Wess and Bagger [26] , except the definition of the Dirac spinor 
as in Equation ([5]). Many useful formulas can be found in Appendix B of the reference. 
The propagators of the left- and right-handed quarks are respectively given as 



.u 



and 



(56) 



{px - fJ'Sxfif 



{px - /^5a,o)^ 



We define the charge conjugation matrix 




(57) 



which satisfies C-^^C = -Y^ and = 



C. 



24 



The quark and anti-quark on-shell projectors are defined as 



1 „x ^ . . 1 



Pl±{p) := 2 (1 ± P^^'""") ' •= 2 ^ ^^■'''^ ) ' (^^^ 

where p := p/lp]. Notice that Pl± = Pr=p and Pr-i- = Pl^- 

We list the useful formulas for inverting the inverse of the Nambu-Gor'kov propagator. 
For n X n matrices A, B, C and D, we have 



'A B\ ^ _ fA~^ +A~^BS^^CA-^ -A-^BSJ^ 
Dj -{ -S^'CA-^ S^' 

provided that the matrices A and Sa ■= D — CA^^B are invertible. We also have 



A B\ ^_(-C-^DSc^ C-^ + C-^DSc^AC-^\ 
CD) -\ -Sc'AC-' ) ' ^^"^ 

provided that C and 5c := B — AC'^D are invertible. Other convenient formulas are the 
following. 

{Pu - f^6^,oK = {po + \p\- ^J)(J^PR+ + {po -\p\- fJ')(r'^PR- , 

{Pu - /x(^.,o)^" = {Po + \P\- f^)^'PL+ + {Po -\P\- /x)aOPL_ . (61) 

{[p,y - fidiyfija I = h 



Po-\p\- P Po + \p\- fJ'' 
u X (^°Pr+ , ct^Pr- 

{{p^ - fid^^fija } = h —= . (62 

Po-\p\- P Po + \p\- fJ' 

Pl± o'^Pr± = , Pl± ct'^Pr^ = cj'^Pr^ , 

Pr± a''PL± = , Pr± a'^PL^ = u^Pl^ . (63) 

B Solving Gap Equation 

We solve Equation ([52]) . following References [28l[19]. We split the integration range into, 
< Po < ^0 s-iid /cq < Po < 2fi. Then the dominant contributions in the logarithms of the 
integrand are In ^r?^ = 2 In — In — for the former integration region and In — for the 
latter. We introduce the parameters 

x:=ln^, y:=ln— , xq := In , (64) 
ko Po ^0 

where Sq := S(0). Then the gap equation is cast into the form 

nx) = f^^^\2xj ny)dy- j yny)dy + j^yny)dyy (65) 
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Note that we have S(xo) = — S(0). We can take the derivative with respect to x twice to 
get 

S" = - 



47r2 Nr 



and the appropriate solution is 



S(x) = Socos W^^^J , with ^0 = ^^^^^ 



From the definition of xq, we obtain 

So = exp 



'47r2 AT^ 

^2 _ 1 



(66) 



(67) 



(68) 
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